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We employ a simple multiple scattering model to inves- 
tigate the inclusive reaction p + A p + X for projectile 
momenta in the 100-200 GeV/c range. We find that data 
are consistent with a class of interaction models in which the 
stopping power of nuclei is rather low. We discuss extrapola- 
tion to ultra-high energy and the application to interpretation 
of cosmic-ray air showers at energies up to 10^" eV. 
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I. INTRODUCTION 

To explore the cosmic-ray spectrum beyond about 
10"'^^ eV requires ground-based experiments with large 
effective area and long exposure times to overcome the 
increasingly low flux implied by the steeply falling energy 
spectrum. Such air shower experiments cannot observe 
the primary particle directly but can only sample the 
cascade it generates in the atmosphere. For this reason, 
obtaining results of astrophysical interest, such as the rel- 
ative fraction of different types of nuclei or the fraction of 
gamma-rays, requires extensive Monte Carlo simulation 
to model the cascades and interpret the measurements. 

A recurring problem is that uncertainties in the input 
to the calculations introduce corresponding ambiguities 
in the interpretations of the experiments. A major, and 
to some extent unavoidable, source of uncertainty is the 
modelling of the hadronic interactions at energies well be- 
yond those explored at accelerators. In the extreme case, 
protons with energies near the Greisen-Zatsepin-Kuz'min 
(GZK) cutoff § at ~ 5 X 10"'^^ eV correspond to center-of 
mass energies more than two orders of magnitude beyond 
the highest energy hitherto available at hadron colliders. 
There are uncertainties also in the region of the "knee" of 
the spectrum around 5 x 10"'^^ eV even though this is ap- 
proximately equivalent to the center of mass energy of the 
Fermilab collider. There are two reasons that significant 
uncertainties remain even at this relatively low energy. 
One is that interactions in the atmosphere involve nuclear 
targets, and in some cases nuclear projectiles as well. The 
other is that it is the forward fragmentation region of the 
collision — largely unexplored by collider detectors — that 
primarily determines the rate of energy deposition that 
generates the core of the atmospheric cascades. 

The most global properties of minimum-bias hadronic 
interactions determine the development of air showers. 



These include the cross section and the inelasticity. The 
pp cross section is directly measured up to -^/s w 2 TeV, 
and its extension to higher energy can be obtained by 
extrapolation of fits based on Rcgge theory. |^ There 
is less agreement on how to extrapolate inelasticity and 
related quantities that determine the rate at which en- 
ergy is deposited in the atmosphere via electromagnetic 
subshowers. We focus on inelasticity in this paper. 

We are motivated to study this problem now because of 
intense experimental activity and ambitious new propos- 
als aimed at the highest energy cosmic rays |-§ as well 
as highly instrumented hybrid arrays aimed at discover- 
ing the sources of cosmic-rays that give rise to the knee 
feature in the spectrum, for example |p[-p^. There is a 
corresponding interest in simulations as illustrated by the 
systematic comparison of several codes undertaken by the 
group at Karlsruhe. ||l^] By installing several hadronic 
event generators into the same cascade code, they have 
isolated differences due to the input physics of the inteac- 
tion models from possible technical differences in how the 
cascades are followed. The latter in principle should not 
be sources of uncertainty in any case, being determined 
by well-known physics such as energy-loss by ionization, 
pair-production, bremsstrahlung, etc. 

Qualitatively, the inelasticity of a hadronic interaction 
is the fraction of the beam energy not carried off by the 
fragment of the incoming particle. This fraction of the 
energy is then available for particle production, includ- 
ing neutral pions which transfer energy from the hadronic 
core of the shower into electromagnetic subshowers. In- 
elasticity is just one moment of one of the inclusive distri- 
butions, but it is arguably the most significant for cascade 
development (next to the inelastic cross section itself) be- 
cause it determines the rate at which the initial energy 
of the cascade dissipates. 

For pp collisions there is a precise experimental defini- 
tion that involves the inclusive cross sections for produc- 
tion of protons, neutrons and their antiparticles. Defin- 
ing 
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(because there are two nucleons in the initial state) and 
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Here Kd is the elasticity, and the inelasticity is defined 

as 
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For pion-initiated interactions, a precise definition of 
elasticity requires a model because of the essential am- 
biguity between produced and fragment pions. At high 
energy the elasticity defined in Eq. ^ is approximated 
by the integral over the leading nucleon in the lab frame, 
which is the definition we use in the remainder of this 
paper. 

For p-nucleus collisions we follow the work of Refs. 
and consider partial inelasticities in the frame- 
work of a Glauber multiple scattering formalism [ p^ . 
The elasticity is given by 
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is the probability for exactly v wounded nucleons in a 
target of mass A and (E)v is the mean energy of the 
leading nucleon in collisions with exactly v wounded nu- 
cleons. The partial inelasticity coefficient is defined by 
the relation 
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The total inelastic hadron-nucleus cross section is 
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We calculate both crP^^j and the partial cross sections o-p* 
from the cross sections for pp scattering and the nuclear 
profiles as described in Rcf. |p7|] . 

As emphasized in Ref . , there is no basis for a naive 
interpretation of Eq. because fast fragment (s) of the 
projectile do not reach an asymptotic physical state until 
well outside the nucleus. The strategy is to assume that 
1\ is determined by pp scattering and to treat the remain- 
ing partial inelasticities (y > 1) as free parameters con- 
strained by fitting p-nucleus data within the framework 
of the model. Our method and conclusions are similar 
to those of Ref. , although we have improved on their 
analysis by using a larger data sample, by treating neu- 
trons and protons separately and by considering effects 
of diffraction. 

The paper is organized in the following way. Section 2 
describes in detail the multiple scattering model we use 



for the description of proton-nucleus interactions. Sec- 
tion 3 introduces the fits to pp data that we need for the 
definition of Ii and other parameters for the case i/ = 1. 
In section 4 we give the fits to p-nucleus data and the 
resulting values for /j/>i. Section 5 contains a discus- 
sion of the results in the context of models in current 
use for calculations of cosmic ray cascades at extremely 
high energies, and examples of estimated inelasticities in 
proton air collisions at very high energy. Section 6 gives 
a summary of the results and conclusions. 



II. THE MULTIPLE SCATTERING MODEL 

The outgoing nucleon in the reaction p -|- A ^ N -I- X 
(N being either neutron or proton) can be specified by its 
transverse and longitudinal momenta, pt and x = pi/po 
where Po is the incident proton momentum. We model 
the differential cross section for this process as a sum over 
final state distributions corresponding to definite num- 
bers of wounded nucleons, v, 
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The transverse momentum distributions are assumed to 
be described sufficiently by an exponential form for fixed 
values of x and are specified by the slope functions b^{x). 
Longitudinal momentum distributions for final state nu- 
cleons are contained in the functions M^{x) which are 
normalized as 



and 
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The numbers express the outgoing nucleon multiplici- 
ties for each number of wounded target nucleons. So MP 
and M" give the x-distributions and relative numbers of 
protons and neutrons after v collisions. Eq. ^ expresses 
the fact that our analysis follows only the forward out- 
going nucleon. When Eq. || is integrated over all final 
nucleon momenta one recovers the inelastic pA cross sec- 
tion times the mean nucleon multiplicity as expected. 

Experimental data on the processes p -|- p ^ p -I- X and 
p 4- p — > n -I- X may be used to fix the {x,pt) distribu- 
tions and the nucleon multiplicities for v = 1. For larger 
numbers of wounded nucleons we employ the iterative 
scheme discussed above (Eq. 7), which is similar in spirit 
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to that used by Hwa and also Hufner and Klar. jlf 
The longitudinal distributions are related by 
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The superscripts + and — above describe interactions 
which preserve and change the projectile isospin respec- 
tively, with the parameters (3 specifying the fraction of 
isospin preserving reactions. After v — 1 collisions, a nu- 
cleon having longitudinal momentum fraction x/y has 
probability S^^-^^{y) to transition to a state having mo- 
mentum fraction x. 

For the case of an incident proton beam, we can take 
the nucleon distributions after zero collisions as boundary 
conditions, 



and 
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and require that Eqs. ^ and |Tj reproduce the exper- 
imentally determined distributions Mf{x) and Mf(a;). 
One immediately finds 
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and /3o = n\ 

Generalizing this result to allow for different inelas- 
ticities upon subsequent collisions, we adopt the power 
law form with a set of adjustable parameters, a^, to be 
determined by fits to pA data. 
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When — we recover a 'naive' multiple scattering 
model in which all collisions proceed equally like isolated 
pp events. With > contributions from non-leading 
collisions {i' > 1) are harder than the initial collision. 
We will show that the existing data on pA interactions 
strongly supports leading baryon spectra that are signif- 
icantly harder for v > \ than for v — 1. 

In order to define partial inelasticity within our formal- 
ism we calculate the mean value of x after v collisions. 



< x >P'"= /q^ dxa;MP'"(x) and relate it to < x >l''^^ 
with the ratio giving the elasticity coefficient for the z/"^ 
collision (or one minus the inelasticity coefficient). Inte- 
grating Eqs. |l^ and ^in this manner one arrives at the 
relation, 



[nl <x>l +nl <x>l] = [l-h 



< X > 



(20) 



where the mean inelasticity is 

^ 1 - (3^_i / dyyS^_j^{y) 
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III. FITS TO P + P DATA 

In the current model, the forms determined for 
Mf''^{x) and b^''^{x) by fitting the available p-f p p+X 
and p + p ^ n + X data can be thought of as a set of 
initial conditions that play a crucial role in what we will 
eventually infer about leading baryon inelasticity from 
the p-nucleus data. This is true primarily because the 
1/ = 1 term of Eq. ^ represents 20-30% of the cross sec- 
tion even for the heaviest nuclei, and secondarily due to 
the connection between M^^i{x) and Mf{x) prescribed 
by Eqs. ^ and Our fitting procedure is straightfor- 
ward; we use Eqns with = 1 and expand Mf''^{x), and 
6P'"(a;) each in a finite Taylor series. The coefficients are 
then adjusted to minimize the P^r degree of freedom 
when compared to data. 

Figure (|l|) shows the differential cross section for the 
process p-|-p^p-|-Xasa function of longitudinal mo- 
mentum fraction, x, for several values of the transverse 
momentum between 0.3 and 1.0 GeV/c. The data are 
for beam momenta of 100 and 175 GeV/c taken from ref- 
erences [|l^ and and the curves represent our best 
fit. 

In Fig. (H) we plot the differential cross section for 
p + p^ n + Xasa function of x for four transverse 
momentum bins from 0.15 to 0.75 GeV/c. The neutron 
data correspond to lab system beam momenta of 282, 
500, 1060, and 1500 GeV/c from reference and the 
curves show our best fit. Our fits suggest ^^(a;) ~ b^{x), 
so we will not distinguish between them in the remainder 
of our discussion. 

In fitting the distributions of protons we have sepa- 
rated the single diffractive and non-single diffractive com- 
ponents, so that MP(x) = MP'^'^{x) + MP'^'^^'^ix). The 
forward diffractive component (target dissociation) rep- 
resents approximately 10% of the inelastic pp cross sec- 
tion. For its a;-dependence we use the functional form 
(1 — x)^^ [|l8|, with a kinematical cutoff near x = 1. 
With forward diffraction fixed, wc then fit AfP'"'"^(a;) by 
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FIG. 1. Differential cross section for p + p ^ p + X plotted 
versus momentum fraction x for several transverse momentum 
bins. The data are for beam momenta of 100 and 175 GeV/c 
from references 19) and 20). The curves show our best fit to 
the data. 




the procedure described above. In our analysis of pA 
data below, we make the approximation of including the 
diffractive component only in the case = 1. This means 
that diffractive contributions to pA spectra are confined 
in our approach primarily to the region x > 0.85. 

The functions Aff Mf , and bi obtained from the 
best fits are shown in Fig. (H) along with the diffractive 
component Mf'^*^ to complete the picture. We find for 
non-single diffractive protons < x >^''^^'^— 0.44 and for 
neutrons < x >i~ 0.26. The proton and neutron mul- 
tiplicities derived from our fits (including single diffrac- 
tion) are = 0.62 and = 0.27. This should be com- 
pared with a proton/neutron ratio of 2:1 in a naive con- 
stituent quark picture of non-diffractive collisions. Given 
the overall normalization uncertainties in the pp data 
(estimated to be ~20%), the fitted results are remark- 
ably close. For the calculation of nuclear processes the 
distributions are normalized so that there is exactly one 
nucleon (proton or neutron) propagating through the nu- 
cleus. (For the first wounded nucleon only the normalized 
distribution includes the diffractive component.) 

We note from our fits that final state neutrons appear 
to be significantly softer than their non-diffractive pro- 
ton counterparts, although the uncertainties in the data 
for — > n-f X are larger than for production of pro- 
tons. This difference in momentum distribution for neu- 
trons and protons can have two important consequences. 
First, in Eqs. |l3| and |l^ there is a 'mixing' of neutron 
and proton spectra controlled by the parameters /9iy>o- 
It is easy to see that the amount of this mixing can af- 



FIG. 2. Differential cross section for p -I- p ^ n -I- X plotted 
versus momentum fraction x for four transverse momentum 
values: a) 0.15 GeV/c; b) 0.30 GeV/c; c) 0.50 GeV/c; and d) 
0.75 GeV/c. These ISR data correspond to laboratory system 
beam momenta of 282 (filled circles), 500 (open circles), 1060 
(filled squares), and 1500 GeV/c (open diamonds) from Ref. 
21. We have used fits of the form of Eq. 9 {v=l) to interpolate 
the data to selected values of x and px- The curves show our 
best fit to the data. 
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FIG. 3. Probability distributions for final state protons, 
Mf-'""^{x), Mf''"^{x), and neutrons, Mi{x), extracted from 
fits to p + p data. Also shown is the slope function, b^(x). 
The mean momentum fractions and multiplicities are quoted 
in the text. 

feet fits to p-nucleus proton spectra if the input [ly = 1) 
speetra differ. We will examine this in the next section. 
Second, because the inital neutrons are softer, the lead- 
ing nucleon inelasticity relevant to high energy cascade 
simulation may be less than one predicted on the basis 
of proton data alone together with the assumption that 
the inclusive distribution of neutrons is similar to that of 
non-diffractive protons. 

IV. FITS TO P + NUCLEUS DATA: 
INELASTICITY 

The number of terms to keep in Eq. ^can be guided by 
considering the Glauber probabilities cr^^/afj^i for heavy 
target nuclei (mass ~ 200). One finds that roughly 90% 
of the cross section is obtained by the first five terms and 
99% by the first eleven. We have terminated the sum at 
v = 12. 

Twelve terms results in a large number of parameters 
to be fit unless some additional assumptions are made. 
Our approach is to treat all interactions subsequent to 
the initial one on the same footing. This means that 
Q!y_i, /3^-i, and the functions b,^{x) have the same value 
for > 1. Indeed we have checked that relaxing this 
constraint has no substantial impact on the results. Only 
marginally better fits are obtained if, for example, we 
allow a different a for each value of v. 



FIG. 4. Differential cross section for the process 
p + A — > p + X plotted versus momentum fraction x at fixed 
pt ~ 0.3 GeV/c for 100 GeV/c protons on targets ranging 
from hydrogen to lead. The data are from reference 20). The 
dashed and solid curves show our best fits with 13 — 2/3 and 
(3—1 respectively. 

We give the slope functions the simple quadratic form 
i>v>i(x) = a -\- bx + cx^ and have checked that higher 
order terms do not substantially improve the fits to data. 
These three parameters together with a^^^i and /3^>i give 
a total of five free parameters for fitting the p-nucleus 
data. 

We have examined two distinct cases for the value 
of /3i/>i based on different extreme pictures for the nu- 
cleon propagation through the nucleus. The first can 
be thought of as the naive case in which all interac- 
tions proceede identically; that is, the probability for 
isospin preserving reactions at each step is just equal 
to the proton multiplicity observed in p -|- p reactions, 
fiv>i = '^i — 2/3. Note that in this case the probability 
that the leading nucleon is a proton quickly approaches 
1/2 with increasing v (14/27 for = 3 and is f^^)- The 
opposite extreme is that the isospin of the leading nu- 
cleon is determined solely at the first interaction. This 
second case corresponds to I3v>i = 1. We will show that 
these two pictures lead to somewhat different conclusions 
for the inelasticity of non- leading interactions required to 
fit the data. 

The data we use to study inelasticity are inclusive pro- 
ton spectra from p -I- nucleus reactions by Barton et al. 
pO| and Bailey et al. with beam energies of 100 and 
120 GeV respectively. The 100 GeV data were collected 
for C, Al, Cu, Ag, and Pb targets for two transverse 
momentum bins of 0.3 and 0.5 GeV/c. The 120 GeV 
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FIG. 5. Same as Fig. (4) but with pt = 0.5 GeV/c. 



data were reported summed over transverse momenta for 
Be, Cu, Ag, W, and U targets. These data afong with 
om' best fits are shown in Figs. (||), and (^). The 
dashed and sohd curves in these figures correspond to 
the /3 = 2/3 and /3 = 1 scenarios respectively. 

In Fig. (|^) we show the per degree of freedom statis- 
tic of these fits plotted as a function of the inelasticity of 
non-leading interactions as determined according to Eq. 

The /? = 2/3 and /3 = 1 scenarios yield their best fits 
for I — .14 and .18 respectively. The [3=1 case offers a 
somewhat better overall fit to the data. It is clear that 
a naive multiple scattering picture which corresponds to 
the case of /3 = 2/3 and / = .5 is certainly not supported 
in the present analysis. 



V. IMPLICATIONS FOR CASCADES 

Most hadronic models currently in use for calculation 
of air showers at high energy fall into one of two cate- 
gories. One group [p3|~p6t is based on the dual parton 
model (DPM) |27| or the related Quark-Gluon String 
model (QGS) ]28| . Another approach is to use some 
variation of statistical or thermodynamical ideas |p9| , ^ , 
producing particles via clusters or fireballs, but con- 
strained to agree with the observed persistence of some 
high-energy fragments of the projectiles. There are sev- 
eral examples of this type of model, for example Refs. 
pl|-p3[. Here we focus on one particular model p^ ] 
that has been used recently to reevaluate implications of 
the Fly's Eye measurements for cosmic-ray composition 



FIG. 6. Differential cross section da/dx for the process 
p + A ^ p + X plotted versus momentum fraction x for 120 
GeV/c protons on targets ranging from Beryllium to Tung- 
sten. The data are from reference 22). The dashed and solid 
curves show our best fit with /? = 2/3 and 13 = 1 respectively. 




V > 1 



FIG. 7. Quality if fit, per degree of freedom, is plotted 
versus inelasticity of non-leading collisions, Iv>i, for our fits 
to p + nucleus data. Separate curves for the /3 = 2/3 and 
/3 = 1 cases are shown. The minima occur at Iv>\ = 14% 
and Iv^i — 18%, with the latter providing a somewhat better 
overall fit. 
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TABLE I. Proton-proton and p-air inelastic cross sections 
with corresponding mean number of wounded nucleons from 
Eq. 8. 



Energy (GeV) 


(Tpp (mb) 


(Tp_air (mb) 


W) 


10^ 


33.0 


284. 


1.69 


10^ 


67.0 


427. 


2.29 


10^ 


102. 


542. 


2.77 


10" 


142. 


661. 


3.14 



TABLE II. Inelasticities for proton-proton interactions 
and for proton-air interactions for two classes of models: 
A=string-type models and B=statistical-type models (see 
text). The headings for p-air correspond to different assumed 
values of /^>i, e.g. 7^>i = 0.14, 0.18, etc. 
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B 


h 


.14 


18 


HE) 


.14 


.18 Ji 


10'' 


.50 


.50 
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10^ 


.55 


.26 


.62 
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74 


.58 


.38 


.40 .45 


10^ 


.57 
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.66 
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79 


.60 


.35 


.39 .40 


10" 


.58 


.15 


.68 


.70 


81 


.61 


.35 


.39 .36 



around 10^^ eV. The model of Ref. fsl is an extrapola- 
tion to high energy of the work of Chou and Yang 1 3^ . 

In the first group of models minimum-bias hadronic 
interactions proceed by the exchange of strings stretched 
between fragments of the incoming projectile and target 
particles. Strings radiate a characteristic multiplicity of 
secondaries per interval of rapidity, so the increase of 
multiciplicity is essentially logarithmic in energy (or more 
accurately, a power of the logarithm because the number 
of exchanged strings increases with energy) . Inelasticity 
is determined basically by the momentum distributions of 
the valence constituents, increasing slightly with energy 
as more soft strings (coupled to sea quarks rather than 
valence quarks) are exchanged. 

The cluster models are generally characterized by a 
more rapid, power-law dependence of multiplicity on in- 
variant mass of the produced clusters. The observed rise 
of multiplicity in the central region is then matched by 
requiring the events to become increasingly elastic as en- 
ergy increases, so that the fraction of total event energy 
going into particle production decreases while the frac- 
tion going into the leading nucleons increases. (See Ref. 
p6t for a discussion of inelasticity in the context of this 
class of models.) 

In Table 1 we show estimates for pp and p-air cross 
sections along with the mean number of wounded nucle- 
ons per interaction from Eq. |^ for lab energies ranging 
from ISR to those relevant in EAS analysis. We note here 
that the range of nuclei used in our study of inelasticity, 
Be {A w 9) to Pb {A « 207) and \J {A ^ 238), is nicely 
matched to the range of energies we wish to consider for 
hadron collisions in air. For a nucleus of mass 200, using 
a standard estimate of the p-nucleus cross section, p^ ] 
we find (1^)200 ~ 3.77 from Eq. B at low energy. In com- 



parison, the mean number of wounded nucleons expected 
in a proton-air collision at 10^° eV from Table 1 is 3.14. 

Extrapolations of the two different types of models for 
hadron-hadron interactions beyond collider energies di- 
verge significantly. We illustrate this in the first section 
of Table 2 by listing the inelasticity for pp collisions, /i, 
as a function of energy to represent the two classes of 
models. The pp inelasticity in column A is chosen to be 
similar to that of Ref. |2^ , while that for B is from the 
work of Ding et al. . In both cases we use for illustra- 
tion the traditional value of 0.5 at low energy, rather than 
the somewhat higher value implied by the distributions 
in Fig. 3. 

Next we calculate the corresponding inelasticities for 
p-air collisions in the two classes of models startingfrom 
the the assumed values of Ii and using Eqs. |, |, I and 
1^ to calculate overall inelasticity for various assumption 
about /j/>i. Based on the analysis of this paper, we use 
Ii^>i = 0.14, 0.18. We also show the result of the 'naive' 
model of propagation through the nucleus = h) for 

illustration, although we have seen that it is inconsistent 
with existing data. (The column labelled 1(E) is dis- 
cussed below.) At the highest energies, we see that 'sta- 
tistical' models predict characteristic energy losses per 
collision of only 35 to 40% compared with about 60 to 
70% for 'string' models. Even at energies characteristic 
of the "knee" region, the models are already predicting 
significant differences in energy deposition rates for pro- 
ton initiated air showers. 

The DPM and QGS models incorporate scattering on 
nuclear targets explicitly. For example, when only one 
target nucleon is wounded, a constituent quark(di-quark) 
belonging to the projectile proton couples to a string that 
in turn connects to a di-quark(quark) belonging to the 
wounded nucleon. In cases where there are two or more 
wounded nucleons in the target, the additional nucleons 
are coupled only to the sea quarks of the projectile. In 
this way the desired physics is reproduced by the model. 
In particular, the excited nucleon, being off mass-shell, 
does not interact repeatedly as a physical nucleon inside 
the nucleus. Moreover, the extra multiplicity characteris- 
tic of a collision on a nuclear target is naturally confined 
to the central region and the target fragmentation re- 
gion of phase space. Capella et al. [ p7[ point out that in 
their model the partial inelasticity (X^i) is of order 0.2, 
decreasing slightly on successive collisions in the same 
nucleus. We have checked that SIBYLL also shows 
this behavior. Thus the string-type models are consis- 
tent with the result of our analysis of proton-proton and 
proton-nucleus collisions. In addition, we note that in 
this type of model it may be more natural to make the 
choice P = 1; that is, to assume that the ultimate identity 
of the final state nucleon is determined only once during 
the interaction with the nucleus. 

We note here that the singular nature of the sea quark 
distributions for small momentum fractions leads to a 
threshold effect in string- type models for /i/>i. Asymp- 
totically the sea-quark momentum on the projectile side 
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VI. SUMMARY 




10'' 
V s [GeV] 



FIG. 8. Mean partial inelasticity, I^>i, versus energy from 
the SIBYLL interaction model. The decrease with increasing 
energy is intrinsic to string-type models and is a consequence 
of the singularity in the momentum distribution of sea 
quarks. 



becomes negligible so that the fractional momentum re- 
moved from the projectile by wounded nucleons with 
v > 1 is small at high energy. This leads to a decrease 
in the value of /y>i as energy increases. This behavior is 
characteristic of string- type interaction models. To illus- 
trate, we use the joint probablility distribution for pro- 
jectile partons from the SIBYLL interaction model to 
evaluate the fraction of energy removed from the projec- 
tile for different numbers of wounded nucleons. Fig. (|^) 
shows the result for averaged over different values 

o{ ly > 1. We see that energy losses due to interaction 
with the quark sea of the projectile decrease significantly 
at high energy. The column in Table 2 labeled I{E) uses 
these energy dependent values of /i/>i for the estimate of 
the overall inelasticity. 

The generalization from pp to p-nucleus is not pre- 
scribed in the statistical models, at least not in the form 
used by Ding et al. |Q. Given the observed rapid en- 
ergy deposition in air showers, users of statistical models 
generally adopt the "naive" treatment of inelasticity in 
nuclei {Iu>i = ^i) to compensate for the intrinsically 
high degree of elasticity of the hadron-hadron model. It 
has also been used in the context of some quark models 
of hadron-hadron interactions |^7| , making the hadron- 
nucleus interactions highly inelastic. 



Analysis of particle production in proton-proton and 
proton-nucleus collisions within a multiple scattering 
framework leads to the conclusion that the second and 
higher interactions of the excited nucleon inside the nu- 
cleus are relatively elastic. Assuming this feature of nu- 
clear interactions persists to high energy, we can esti- 
mate the inelasticity in hadron-nucleus collisions beyond 
the energy range for which we have data. The results 
depend on the behavior of the cross section and inelas- 
ticity for proton-proton collisions, as illustrated in Tables 
1 and 2. Since the column labelled Iu>i = Ii is ruled out 
by the p-nucleus data, we conclude that the inelasticity 
on nuclear targets in the statistical models at high en- 
ergy must be quite low. As pointed out in Ref. Q, (see 
also Ref. [^), such low inelasticity is unable to account 
well for the Fly's Eye data. [|0|j39| Models of the type 
QGS and DPM represent interactions on nuclear targets 
in a way that is consistent with the low energy data on 
nuclear targets. They predict that inelasticity increases 
slowly with energy, with a modest increase for nuclear 
targets. 
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